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We examine the adsorption of a single Ni atom on a monolayer of MgO on a Ag substrate 
using DFT and DFT+17 computational approaches. We find that the electronic and magnetic 
properties vary considerably across the three binding sites of the surface. Two of the binding sites 
are competitive in energy, and the preferred site depends on the strength of the on-site Coulomb 
interaction U. These results can be understood in terms of the competition between bonding and 
magnetism for surface adsorbed transition metal atoms. Comparisons are made with a recent 
experimental and theoretical study of Co on MgO/Ag, and implications for scanning tunneling 
microscopy experiments on the Ni system are discussed. 
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I. INTRODUCTION 

The study of magnetic adatoms on surfaces has drawn 
recent attention due to possible applications in the realm 
of magnetic storage and quantum computation. The 
density of magnetic storage has enjoyed exponential 
growth for several decades, but this growth will even¬ 
tually slow as particle sizes approach the superparamag- 
netic regime.^ An alternative, bottom-up approach is to 
start from the atomic limit. With the scanning tunneling 
microscope (STM), single atoms can be moved around 
on a surface to construct desired nanostructures, and the 
STM can also be used to probe the electronic and mag¬ 
netic properties of those nanostructures. Such experi¬ 
ments have found, for example, that a single Fe atom 
on a Cu 2 N monolayer (ML) on Cu has a large magnetic 
anisotropy energy,^ and that a magnetic bit consisting of 
an array of 12 Fe atoms on the same surface has a stable 
moment that can stay in the ‘on’ or ‘off’ state for hours 
at cryogenic temperatures (^ 1 K).^ 

A recent STM study of a single Co atom on a ML of 
MgO on the Ag(OOl) surface found that the Co atom 
maintains its gas phase spin of S = 3/2, and, because of 
the axial properties of the ligand held, it also maintains 
its orbital moment on the surface (L = 3). The resulting 
magnetic anisotropy is the largest possible for a 3d tran¬ 
sition metal, set by the spin-orbit splitting and orbital 
angular momentum. The measured spin relaxation time 
of 200 /xs is three orders of magnitude larger than typical 
for a transition-metal atom on an insulating substrate.^ 

Here we investigate the structural, electronic, and mag¬ 
netic properties of Ni adatoms on the same substrate, 
MgO/Ag(001), using density functional methods. Previ¬ 
ous authors®^^^ have studied the adsorption of transition 
metal atoms on an MgO substrate using ab-initio tech¬ 
niques. These studies, which employed various surface 
models and approximations for the exchange-correlation 
functional, are in general agreement that the preferred 
binding site of a Ni adatom on the MgO(OOl) surface is 


on top of an O atom, and that structural distortion of 
the MgO surface upon Ni adsorption is minimal. The 
situation is less clear when it comes to the spin state of 
the adatom, since the s-d transition energy is so small in 
Ni. Calculations that employ the generalized gradient ap¬ 
proximation for the exchange-correlation functional gen¬ 
erally predict a full quenching of the Ni moment on MgO, 
while partial inclusion of Fock exchange as in the B3LYP 
hybrid functional predicts that the triplet spin state is 
slightly more favorable. 

The system studied in the current work differs from 
previous studies of Ni adsorption on MgO because the 
substrate consists of a single layer of MgO atop the 
Ag(OOl) surface. This aligns more closely with STM ex¬ 
periments that use a thin insulating layer to decouple the 
magnetic adatom from the conducting substrate that is 
necessary for electrically probing the system. We show 
that, in contrast to the surface layer of an MgO sub¬ 
strate, the MgO ML on Ag can deform significantly due 
to interactions with the Ni, resulting in a very different 
potential energy surface for adsorption. Further, we in¬ 
vestigate how the interaction of Ni with the MgO/Ag 
substrate is affected by on-site Coulomb interactions and 
find that the preferred binding site depends on the inter¬ 
action strength U. Unlike the case of Co adatoms on the 
MgO/Ag substrate, the Ni spin moment is always lower 
than that of the isolated atom. The degree to which the 
moment is reduced depends strongly on the binding site. 
These results can be understood by considering the Ni 3(i- 
4s hybrid orbitals that participate in bonding at different 
sites. We conclude with a discussion about experiments 
that could corroborate our findings. 
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II. METHOD 
A. Computational methods 

Two sets of density-functional-theory (DFT) calcula¬ 
tions were carried out, one using the linearized aug¬ 
mented plane wave (LAPW) method as implemented in 
WIEN2K,^^ and the other using the projector augmented 
wave (PAW) method^^ in the VASP package.Both 
started with the Perdew, Burke and Ernzerhof formu¬ 
lation of the generalized gradient approximation (GGA) 
for the exchange-correlation functional.^® In the LAPW 
calculations, structures for each binding site were relaxed 
within the GGA. These structures were then used to 
calculate Hubbard U values for each site using the con¬ 
strained DPT method of Madsen and Novak.GGA-I-G 
calculations were then carried out using those site-specific 
values of U to examine electronic and magnetic proper¬ 
ties, without further structural relaxation. On the other 
hand, PAW calculations were used to optimize structures 
within both GGA and GGA-|-[/. However, since the total 
energy in DET-I-G methods depends on the value of the 
on-site Goulomb interaction strength, the same value of 
U was used for each adsorption site to allow comparison 
of binding energies. In all GGA-I-H calculations, a ro- 
tationally invariant method in which only the difference 
U — J is meaningful was employed.^® 

Within GGA, the two sets of calculations yielded very 
similar results for adsorption geometries, binding ener¬ 
gies, and electronic and magnetic structure, as expected. 
Although the calculations were used for different pur¬ 
poses in assessing the effect of on-site Coulomb repulsion, 
the GGA-I-G results from the two methods were generally 
consistent and displayed the same trends.^® 

B. Supercell geometry and binding sites 

STM experiments require a conducting substrate, yet 
for magnetic nanostructures, it is desirable to suppress 
interaction between the adatoms and a metallic sub¬ 
strate. Hence the adatom is often placed on a thin 
insulating layer at the surface of the substrate. It has 
been shown experimentally that ultrathin ionic insulat¬ 
ing layers can shield the adatom from interaction with the 
surface.^® In a recent experimental study of Co adatoms, 
a single atomic layer of MgO was used as the insulating 
layer above an Ag substrate.The lattice constants 
of MgO and Ag are well matched (4.19 A and 4.09 A, 
respectively), and DPT calculations of a single layer of 
MgO on the Ag(lOO) surface have found that is it ener¬ 
getically favorable for the O atoms in the MgO layer to 
sit above the Ag atoms."^ 

Here we used this alignment of MgO on Ag(OOI) in 
inversion symmetric (001) slabs of at least five Ag lay¬ 
ers sandwiched between MLs of MgO. We found only 
minor differences in results for slabs containing five ver¬ 
sus seven Ag layers. The in-plane lattice constant was 



(a) O site (b) Mg site (c) hollow site 


FIG. 1. (Color online) Top view of the three binding sites. 
Cyan spheres represent Mg, small red spheres O, and the yel¬ 
low spheres Ni. 

fixed at the bulk Ag value, and an isolated Ni adatom on 
the MgO/Ag surface was modeled using a 3/^/2 x 3/-\/2 
supercell of the slab with one Ni atom on each surface, 
corresponding to a lateral separation of 8.68 A between 
adatoms. In the out-of-plane direction, the supercells 
contained 7 to 8 layers of vacuum. 

Three high-symmetry binding sites on the MgO surface 
were considered, as shown in Pig. I: Ni on top of an O 
atom, Ni on top of an Mg atom, and Ni above the center 
of the square formed by nearest-neighbor Mg and O sites. 
These will be referred to as the O site, the Mg site, and 
the hollow site, respectively. 

III. RESULTS & DISCUSSION 
A. Binding energetics and geometries 

Since GGA and GGA-hU do not fully describe impor¬ 
tant correlation effects in the isolated Ni atom, the calcu¬ 
lated binding energies are not as reliable as differences in 
binding energy between different adsorption sites. Pig- 
ure 2 shows the total energy relative to that of the O 
site. Within GGA, the O site is slightly favored over the 
hollow site, while the Mg site is about I eV higher in 
energy. Por comparison, we also plot our results for Ni 
on an MgO substrate. Similar to previous reports,we 
find that on the MgO substrate, Ni clearly favors the O 
site, with the hollow and Mg sites lying about I and 1.7 
eV higher in energy, respectively. 

The difference between the potential energy surfaces 
for Ni adsorption on MgO/Ag versus on MgO can be at¬ 
tributed to the greater freedom that the MgO ML has 
to deform to accommodate the adatom. The pure MgO 
substrate remains very flat upon Ni adsorption. When 
Ni is on the O or Mg site, the displacement of surface 
atoms is negligible, and when it is on the hollow site, 
the neighboring O atoms displace out of the plane by 
about 0.09 A. The situation is very different for Ni on 
MgO/Ag. Table I lists the distance between the adatom 
and its nearest Mg and O neighbors, as well as the verti¬ 
cal displacement A^, of those Mg and O atoms. For each 
adsorption site, at least one of the neighboring atoms is 
displaced vertically by more than O.I A. The deformation 
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FIG. 2. (Color online) Total energy per Ni adatom on 
MgO/Ag, calculated for structures relaxed within GGA and 
GGA+17 ([7 = 4 eV) from PAW calculations. For comparison, 
energies calculated for Ni adatoms on MgO are also shown. 
All energies are plotted relative to the O site energy. 


of the MgO ML is most striking when Ni is on the hollow 
site, with neighboring O atoms being pulled up out of 
the MgO ML by nearly 0.5 A. This allows the hollow site 
to be competitive in energy with the O site on MgO/Ag, 
in contrast to what happens on the MgO surface. 

To examine the influence of local correlations on 
the binding energetics of Ni on MgO/Ag, we present 
GGA+U results. Since total energies can only be com¬ 
pared for calculations that use the same value of t/, 
structures were relaxed assuming U = A eV for all sites. 
GGA-hC/ predicts that the O and hollow sites reverse or¬ 
der, as shown in Fig. 2. The site above an Mg atom 
still remains much less favorable than either the O or the 
hollow site. Of course the Coulomb interaction U de¬ 
pends on the local environment of the Ni atom, so it is 
an approximation to assume a fixed value for all adsorp¬ 
tion sites. Nevertheless, these results demonstrate that 
on-site correlations can change the preferred adsorption 
geometry. A similar effect has been reported for the ad¬ 
sorption of Sd transition metals on graphene. 

For a Ni adatom above an O atom, the Ni-0 bond 
is weakened by on-site correlations, as evidenced by an 
increase of about 8 % in the Ni-0 bond length in going 
from U = 0 to C/ = 4 eV. On the hollow and Mg sites, 
differences between geometries relaxed within GGA and 
GGA-bC/ are much smaller. 


Given the limitations of the DFT-I-C/ method, we are 
not able to definitively predict whether the O site or hol¬ 
low site is preferred. However, it is likely that the Mg site 
is uncompetitive, and not experimentally feasible. Hence 
the remainder of the paper deals primarily with the hol¬ 
low and O binding sites, comparing their electronic and 
magnetic properties. 

B. Electronic and Magnetic Properties 

Based on the geometries relaxed within GGA, our con¬ 
strained DFT calculations yield Coulomb parameters U 
= 4.6, 6.0, and 5.0 eV for Ni Sd electrons on the O, hol¬ 
low, and Mg sites, respectively. The GGA-I-C/ results 
presented in this section use these site-specific values of 
[7, but assume the GGA-relaxed geometries. As men¬ 
tioned, the 0 -site geometry is somewhat sensitive to the 
inclusion of on-site Coulomb repulsion, but the Mg-site 
and hollow site geometries are not. 

Atomic Ni has a ground state configuration of 3d®4s^ 
(^^ 4 ) with the next highest state 3d®4s^, belonging to a 
different multiplet i^D), only 0.025 eV away.^^ The av¬ 
erage energies of these multiplets, however, have 3d®4s^ 
0.030 eV lower in energy.^^ The proximity of these atomic 
states sets the stage for a competition between chemical 
bonding and magnetism.® 

Our GGA calculations give a 3d®4s^ ground-state con¬ 
figuration for atomic Ni, which underrepresents the s oc¬ 
cupation compared to experiment, while inclusion of a 
local Coulomb interaction {U = 6.0 eV) yields a 3d®4s^ 
configuration. In both cases, the Ni moment is calculated 
to be 2.0 fiB- 

The calculated magnetic moments for Ni on MgO/Ag 
are listed in Table 11. The Ni spin moment is reduced 
from the atomic value on all sites, though by varying 
amounts. When Ni is on the O site, the magnetic mo¬ 
ment is calculated to be zero within GGA, but it in¬ 
creases slightly upon inclusion of U. (When the structure 
is relaxed within GGA-I-C/ and the Ni-0 bond length in¬ 
creases by 8 %, the moment increases to about 0.3 fiB-) 
On the hollow site, the magnetic moment of about 1 ^b 
is insensitive to U, while the larger moment of about 1.5 
/is on the Mg site increases with U. 

In Table HI, we present the results of a Bader charge 
analysis^^ (space filling) for the adatom and its nearest 


Q .Q Ni on MgO/Ag, GGA 
A A Ni on MgO/Ag, GGA+U 
• ' Ni on MgO, GGA 



TABLE I. Relaxed distance between adatom and nearest- 
neighbor Mg and O atoms, and height of these atoms above 
the MgO ML. All distances are in A. Results were obtained 
within the GGA using the all-electron LAPW method. 


TABLE 11. Total magnetization per adatom (in ps) for dif¬ 
ferent sites, calculated within GGA and GGA-I-C7, from all¬ 
electron LAPW calculations. Site-specific values of U are 
listed in eV. 


Adatom 

Binding site 

dad-O 

<^ad-Mg 

Azo 

AzMg 

Adatom 

Binding site 

m-GGA 

maaA+u 

U 

Ni 

0 site 

1.79 

2.87 

-t0.15 

-t0.06 

Ni 

0 site 

0.00 

0.16 

4.6 

Ni 

hollow site 

1.90 

2.53 

-tO.48 

-0.18 

Ni 

hollow site 

1.04 

1.06 

6.0 

Ni 

Mg site 

3.63 

2.73 

-0.15 

-t0.20 

Ni 

Mg site 

1.40 

1.66 

5.0 

Go 

0 site 

1.85 

2.95 

-0.06 

-0.20 

Co 

0 site 

2.68 

2.78 

6.9 


















4 


TABLE III. Net charge of Ni atom and its nearest-neighbor 
O and Mg atoms according to Bader analysis. Ag/MgO is 
the bare slab without adatoms. These results were obtained 
from all-electron LAPW calculations. 



Ni 

GGA 

Onn 

Mg„„ 

Ni 

GGK+U 

Onn 


O site 

-0.17 

-1.50 

-tl.72 

-0.13 

-1.52 

-tl.71 

hollow site 

-tO.38 

-1.51 

-H.71 

-tO.43 

-1.54 

-tl.71 

Mg site 

-0.26 

-1.62 

-tl.68 

-0.40 

-1.61 

+1.69 

MgO/Ag 


-1.64 

-tl.71 





neighbor Mg and O atoms, as well as for a bare slab of 
MgO/Ag. Because this analysis includes the interstitial 
electrons, these results can provide insight on the nature 
of the bonding between the adatom and the surface. On 
the O site, the Ni adatom gains electrons from its O 
neighbor (which is less negatively charged than it would 
be on bare MgO/Ag). In contrast, on the hollow site, 
Ni loses electrons to the underlying Ag substrate. These 
charge transfer results are relatively insensitive to local 
Coulomb interactions. 

To understand the site dependence of the Ni moment, 
we examine the bonding and electronic density of states, 
starting with Ni on the O site within GGA. While the 
transfer of electrons from the substrate to a 3d®4s^ con¬ 
figuration can only reduce the spin moment, the gain of 
^ 0.2 is not enough to account for the full quenching 
of the moment. The fact that most of the electrons are 
transferred from the neighboring O atom, and the fact 
that Ni adsorption causes the O atom to displace upwards 
out of the MgO ML (Table I), suggest the formation of 
a covalent bond. Indeed, previous studies of Ni on MgO 
discuss a bonding mechanism involving Ni hybrid orbitals 
of 4s and character. of the two orthogonal 
orbitals, which we denote as As+dz ^, is oriented perpen¬ 
dicular to the surface and interacts strongly with O Pz 
orbitals, forming a bonding combination (mostly O Pz) 
about 6 eV below the Fermi level {Ep) and an antibond¬ 
ing combination (mostly Ni As+dz^) about 1 eV above 
the Fermi level. This can be seen in the electronic density 
of states plotted in Fig. 3(a). The other hybrid orbital, 
denoted As—dz^, is concentrated in the plane parallel to 
the surface, and hence interacts much more weakly with 
the surface. This planar hybrid orbital lies just below the 
other Ni d states, about 1 eV below the Fermi level. The 
covalent interaction between Ni and O pushes the anti¬ 
bonding orbital high enough to be unoccupied in both 
spin channels, so the other d and s-d hybrid orbitals are 
fully occupied, yielding a zero net spin. 

When the orbital-dependent Goulomb interaction U is 
included within GGA-I-G, the density of state changes 
significantly since the d and s states are affected differ¬ 
ently. The hybridization between Ni dz^ and s orbitals 
weakens, and the interaction between these orbitals and 
O p states is affected. The partial density of states plots 
in Fig. 3(b) show the majority and minority spin dz 2 


orbitals split by about 4 eV (^ U). The antibonding Ni 
- O orbital, now mostly composed of Ni s states, has a 
small spin splitting of the opposite polarity: the spin- 
down states lie just below Ep and the spin-up ones lie 
just above Ep. Hence even though U induces a moment 
in the d shell within the Ni muffin tin, this moment is 
screened by the oppositely polarized antibonding orbital 
formed from Ni 4s and O 2pz orbitals. This is illustrated 
in Fig. 4(a), which shows a majority-spin (red) isosurface 
with dz 2 character on the Ni atom, surrounded by the 
antibonding Ni s - O minority-spin (blue) isosurface 
(corresponding to minority spin states near ~ —0.5 eV 
in Fig. 3(b)). The total magnetization of the system 
remains close to zero. 

The situation differs for the hollow site. The presence 
of the Ni at the hollow site disrupts the MgO bond net¬ 
work, and the two O atoms that neighbor the Ni atom 
are displaced out of the surface by nearly 0.5 A to form 
bonds with the Ni, as is evident in Fig. 4(b). The re¬ 
sulting polar covalent Ni-0 bond length is only about 6% 
larger than the corresponding 0-site bond (Table I). The 
hollow site has C 2 v symmetry, meaning that the dxz,dyz 
degeneracy is broken. The lobes of the dyz orbital point 
towards the O neighbors and interact with O Py orbitals. 
However, since dyz transforms differently than 4s, sym¬ 
metry considerations prevent the dyz orbital from form¬ 
ing diffuse hybrids with 4s analogous to the AsEdz^ or¬ 
bitals on the O site that enhance the orbital overlap be¬ 
tween the adatom and surface atoms. While symmetry 
permits hybridization between Ni s and dx 2 _y 2 orbitals, 
the hybrid that interacts strongly with O p states has 
mostly Ni s character, while the orthogonal hybrid or¬ 
bital is mostly Ni dx^-y 2 . Hence the Ni d orbitals all 
retain a higher degree of localization than they do on the 
O site. Within GGA, it becomes favorable for the anti¬ 
bonding Ni dyz -Op orbitals to become spin polarized 
and this largely accounts for the total moment of about 
1 pb- This can be seen in Fig. 3(c). 

Because the effect of U in the hollow site case is to 
increase the splitting between the dyz majority and mi¬ 
nority spin states, the occupation of s and d states and 
the total moment are not significantly impacted, as can 
be seen in Fig. 3(d). Figure 4(b) shows the spin density 
for this case. The majority-spin isosurface in red has an¬ 
tibonding Ni dyz - O Py character. While it is hard to 
see from the angle shown, the minority-spin isosurface in 
blue has the character of a As—dx' 2 -y -2 hybrid. 

We return now to the question of why the local 
Goulomb interaction stabilizes the hollow site relative 
to the O site. When Ni adsorbs to the O site it forms 
a strong covalent bond with O through diffuse sd hy¬ 
brids that overlap significantly with O p orbitals. Since 
the orbital-dependent U disrupts the sd hybridization, it 
weakens the binding of the Ni to the O. On the hollow 
site, geometric and symmetry considerations prevent the 
formation of diffuse sd hybrids that overlap strongly with 
O, so the impact of U, both in terms of binding energy 
and spin moment, is smaller. As a result, increasing the 







5 


GGA GGA+U 



-6 -4 -2 0 2 -6 -4 -2 0 2 

E-Ep (eV) E-Ep (eV) 


FIG. 3. (Color online) Densities of states from inside the muffin-tins of all-electron LAPW calculations, (a) GGA results for 
Ni on O site, (b) GGA-|-f7 results for Ni on O site, (c) GGA results for Ni on hollow site, (d) GGA+U results for Ni on hollow 
site. In all graphs, majority spin is plotted on the positive ordinate and minority spin on the negative ordinate. Total DOS has 
been scaled down significantly, and includes contributions from the interstitial region. The relative scales of nearest-neighbor 
Mg and O, as well as Ni s and d DOS have been adjusted to enhance viewability. 


strength of the local Coulomb interaction decreases the 
stability of the O site relative to the hollow site. 

The presence of magnetization in adatoms depends 
largely on the relative strengths of the covalent inter¬ 
action with surface atoms and the intra-atomic exchange 
coupling for the adatom. In a previous study, ^ our cal¬ 
culations found that Co adsorbs to the O site of the 
MgO/Ag surface. We also found that, in contrast to Ni, 
it maintains a magnetization close to the atomic value 
(Table II). This can be explained by differences in the Ni- 
O and Co-0 interactions. Unlike Ni, Co remains charge 
neutral when adsorbed, according to a Bader analysis, 
and the neighboring O and Mg atoms displace down¬ 
ward rather than being attracted to the adatom (Table 
I). Since the s-d transition energy is larger in Co than 


Ni, the hybridization between s and d orbitals is reduced, 
leading to less interatomic overlap and weaker covalent 
interactions. The d orbitals remain localized enough to 
support magnetism. 


IV. CONCLUSIONS 

Using first-principles calculations, we have found a 
strong site dependence of the electronic and magnetic 
properties of Ni on MgO/Ag. In particular, the magnetic 
moment varies from S = 0 to nearly 1 as Ni is moved away 
from surface O. This is in stark contrast with the situ¬ 
ation of Co on MgO/Ag, in which the O site facilitates 
the preservation of Co spin and orbital moments on the 


































































(a) O site 


(b) hollow site 


FIG. 4. (Color online) Spin density isosurfaces of Ni adsorbed on MgO/Ag, from GGA+17 all-electron LAPW calculations. 
Red and blue surfaces correspond to majority and minority spin densities of 0.00175 e/a.u.^, respectively, (a) O site: the 
majority spin d^2 is screened by the antibonding Ni s - O p orbital, (b) hollow site: the antibonding Ni dyz -Op orbital is spin 
polarized. The y axis is parallel to the line joining the two O atoms that neighbor the Ni. These two O atoms are vertically 
displaced out of the MgO plane. O atoms are represented by small red spheres. 


surface. 

Furthermore, we see that the energy surface for ad¬ 
sorption of Ni on MgO/Ag is very different from that 
for Ni on MgO. On MgO/Ag, the O site and hollow site 
are competitive in energy, and it is not clear from theory 
alone which site will be realized in experiments. How¬ 
ever, the properties of the Ni adatom on the two sites 
differ significantly, not just in the magnetic moment, but 
also in the charge transfer and bonding. 

The energetic proximity of the hollow and O sites raises 
the interesting possibility of being able to select one site 
or the other, and tune the magnetic moment by modifi¬ 
cation of the substrate through, for example, application 
of strain or introduction of defects. It would also be in¬ 


teresting to compare bilayer and monolayer MgO on Ag 
experimentally. The bilayer would be much more rigid 
than the monolayer, and hence we would expect the O 
site to be clearly favored for Ni adsorption. 
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using a plane-wave cutoff of 500 eV, k-point sampling of 
8x8x1 grids, and a Gaussian smearing width of 0.02 eV. 
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